Local distortions in slow moving vortex lattices 
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We observed slow moving vortex lattices in NbSe2 using a home built STM at a temperature of 
4.2 K and magnetic fields in the range of 250 — 500 mT. The vortex lattices move coherently at speeds 
of 1 — 10 pm/s due to a slow decay of our magnetic field of 0.5 mT/day. We observe collective speed 
variations which indicate pinning/depinning events. Furthermore, we found - to our knowledge for 
the first time - small and local deviations from the coherent behavior. Most noticeable were local 
lattice distortions with displacements of the vortices of ~ 2 nm from their ideal lattice position. We 
compared the observations to a 2-D simulation of a vortex lattice moving under similar conditions 
and found similar distortion effects near pinning sites. Thus, studying the distortion patterns could 
lead to a better understanding of the vortex-defect interaction. 



I. INTRODUCTION 

NbSe2 is a widely studied material. It is regularly used 
as a test sample for low temperature scanning tunneling 
microscopes operating in magnetic fields since it is easy 
to prepare and shows a charge density wave. In addition 
it is a type II superconductor allowing the calibration of 
the magnetic field by observing the vortex pattern. At 
low magnetic fields pinning centers, such as structural 
defects in the NbSe2 crystal get populated. Towards 
higher fields the Abrikosov latticei is formed. For our 
conditions of field and temperature the lattice should be 
in the Bragg glass phase^i^ with local hexagonal but no 
long range order. It has also been previously reported 
that the vortices move under the influence of a changing 
magnetic field. The observations were performed using 
Bitter technique^ at low fields, using muon scattering^, 
as well as using STM at higher fields^. However, the 
measured velocities were typically several orders of mag- 
nitude higher than in our study (/zm/s vs. pm/s), making 
local details of the vortex motion harder to observe. In 
our case, the vortices move about 0.25 — 2.5 nm between 
two passes of the STM tip resulting in a highly resolved 
time series. This enables us to determine details of the 
vortex motion unseen in previous experiments. 

As described elsewhere^, consecutive images of the vor- 
tex lattice were taken by measuring variations in the local 
density of states near the superconducting gap (typically: 
14ias = 3 mV, It =0.1 nA). The images were automati- 
cally processed to extract the vortex positions with sub- 
pixel resolution (See Fig.[T]). The position of each vortex 
was followed through the image series to extract the ve- 
locity and the vortex lattice was reconstructed to find 
the local lattice constant. We found collective variations 
in the vortex velocity which primarily demonstrate the 
cohesive behavior of the vortex lattice and whose source 
lies thus lies outside of the observed area. In addition 
we found local lattice distortions on the order of ^ 2 nm 
which are the main topic of this paper. Detailed analy- 
sis of the data revealed that the lattice distortion is also 
accompanied by local variations in the vortex velocity. 




FIG. 1: (Color online) Normalized dl/dV curves of the super- 
conducting phase and the vortex core (left) and dl/dV image 
of the vortex phase (400 nm x 400 nm) in various stages of 
being processed (right): a) original, b) inverted and filtered, 
c) threshold to find vortex, d) product of b) and c) to find 
vortex center as a weighted sum. 



which is necessary to maintain the vortex lattice. While 
clearly visible, these higher order effects suffer a poor 
signal to noise ratio and shall not be discussed in this 
paper. The distortions are distinctly different from dislo- 
cation (eg. recent overview of Bitter decoration^) in that 
the lattice ordering and six fold coordination stays intact. 
A second distinction is that we are observing a moving 
vortex matter. Hence, the stationary distortion must be 
rooted in interactions with the host material. In simula- 
tions similar distortion patterns were found near the loca- 
tion of weak pinning centers. While for a detailed study 
of vortex-defect interactions Ginzburg-Landau calcula- 
tions similar to works by Maurer et. al^ are necessary, we 
chose a 2D Monte Carlo simulation to cover the dynamics 
of a larger number of vortices. Previous studies using this 
method mostly sought to explain transport phenomena 
by examining vortex-defect interactions of moving vortex 
matter-i2r22. Another study examined static, metastable 
vortex configurations around a defect site^^. In our case 
we are interested in the dynamic behavior and in com- 
paring the simulation to our data. 

We employed two different methods for measuring the 
distortions as described in section II. Section HI discusses 
the results of our 2D simulations showing similar distor- 
tions. This gives us confidence that the observations de- 
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scribed in section IV arc not experimental artifacts. 



II. DISTORTION DETECTION 

Observing spacial distortions using an STM has to 
be conducted carefully to ensure that the observations 
are not just due to piezo artifacts. In our case, we 
measure lateral distortions on the order of 2 nm us- 
ing STM images with a pixel resolution of typically 
(400 nm)/(128 pixel) = 3.125 nm/pixel. The computer 
algorithm to detect the vortex position uses on average 
77 pixels per vortex. That reduces the uncertainty by a 
factor of 8.8 to 0.36 nm. 

We used two methods to measure the distortions 
within the vortex lattice. In both cases, we first de- 
termined the spacial relation for each vortex within one 
frame. Assuming a triangular lattice, we found the six 
nearest neighbors for each vortex noting the vortex- 
vortex distance. The lattice was then pruned by exclud- 
ing vortices as possible neighbors at distances larger than 
1.5 times the average lattice constant. This leaves some 
vortices, especially close to the edge of the scan area, 
with coordinations less than six. The local layout is used 
for the subsequent calculations. 

Fig. [2] shows both methods applied to the same frame. 
The first method calculates the lattice distortion as the 
displacement of a vortex with respect to the center of 
mass of the hexagon formed by its six nearest neighbors. 
The center of mass (cx, Cy) is calculated using a sum for- 
mula for an N-sided polygon: 
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FIG. 2: (Color online) Measuring the vortex lattice distortion 
by means of center of mass (Method 1) and using an ideal 
lattice (Method 2). Method 1 gives only results for six fold 
coordinated vortices (green circles) whereas Method 2 allows 
to measure the distortion up to the edge of the scan area. 
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FIG. 3: Method 1 of the distortion detection. Six neighbours 
of a vortex form a hexagon (1-6). The center of mass (cx, Cy) 
is used as a reference to determine the displacement (arrow). 



a fixed lattice constant and orientation yielding a set of 
expected coordinates {(ex, Cy)} per image. The measured 
lattice positions {(px,Py)} are mapped onto the expected 
coordinates using a 2-D, 2nd order polynomial for the x- 
and y-coordinate, respectively. The set of parameters of 
the polynomial {(aj,k,hj^k)} are estimated by solving: 



with the area 
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(xi,yi) denotes the zth point of the polygon (i.e., mea- 
sured vortex positions) . Fig. [3] shows a schematic dia- 
gram of the process. This method is easy to implement, 
but it has two drawbacks: First, if one vortex is displaced, 
it influences the estimated position for its six neighbors 
by ^ of its displacement. Second, each vortex has to 
have six neighbors to be considered (circles in Fig. [2]) 
Consequently, vortices close to the edge of the scan area 
can not be evaluated (crosses in Fig. [2]). Nevertheless, 
this straight-forward method does not require any further 
data processing and allows us to validate results obtained 
by other methods. 

The second method is more elaborate. A backtracking 
algorithm retraces the vortex lattice of a given frame with 



The parameters are then used to map the measured coor- 
dinates onto corrected ones {(cx,Cy)}, effectively remov- 
ing distortions introduced by the scan piezo. The dif- 
ference between the corrected and expected coordinates 
shows the lattice distortion: d = c — e. This method 
allows us to consider more vortices and suffers less from 
back-action of the displacement of some of the vortices. 

We find a good agreement when comparing the two 
methods of determining the lattice distortion for three 
series of measurement. The correlation coefiicients are 
summarized in Table ID We find a very good correlation 
for the X- and y-distortion with values well above 0.8. 
Thus, we focus in the following discussion solely on the 
data determined by the second method. 
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Series 


iV 


X 


y 


angle 


Ml 


1 


47091 


0.86 


0.87 


0.76 


0.75 


2 


20349 


0.88 


0.87 


0.77 


0.74 


3 


44436 


0.87 


0.93 


0.87 


0.80 



TABLE I: Correlation of the lattice distortion calculated with 
two methods within each data set. A'': number of points used, 
a;,i/,angle,|d|: correlation of x- and y-component, direction 
and magnitude of distortion, respectively 



Interaction 


Force constant [nN/m] 


Radius [nm] 


Global force 


0.01 - 20 




vortex- vortex 


30 - 40 


200 


vortex-defect 


-10 --1 


20 - 120 



TABLE IL Typical force and distance constants used in the 
simulation. 
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FIG. 5: (Golor online) Distance dependence of the radial com- 
ponent of the lattice distortion and the magnitude of the ve- 
locity next to a defect. The plot the average over a 1.2 nm 
wide ring and time at a given distance as well as the up- 
per/lower bound. An exponentially decaying function was fit 
to the lattice distortion. 



III. SIMULATION 

The simulation of a moving vortex lattice will be dis- 
cussed elsewhere in greater detaiU^. Here it primarily 
serves as a test whether the effect observed in the data is 
at all feasible in order to rule out strange behavior of the 
scan piczo, the microscope electronics, or an influence of 
local topographic features on determining the vortex po- 
sition. Similar to previous work of others, we calculate 
the flow of an ensemble of vortices over a static landscape 
of defects in two dimensions with periodic boundary con- 
ditions under a constant driving force. The full vortex- 
vortex interaction is taken into account. The equation of 
motion for the «th vortex is given by: 

■qVi = '^Fy {X^ - Xj) + Fd (fj - foj) + -^global 
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FIG. 4: (Color online) Lattice distortions (left) and velocity 
(right) distribution in the vicinity of a single defect. The 
arrows represent averages of a 5 x 5 nm^ area of about 3 data 
points each. The scaling factors for the arrows are md = 4 
and my = 0.6 ^^s ' respectively. 



Here 

^V(r) = /v,o-^-^i(|r|/rv.o) 
\r\ 

and 

FD(r) = /D.o- — •e-(l'^/''°-°)' 

7'D,0 

are the vortex- vortesi^ and vortex-defect interactions, 
respectively. Ki is the modified Bessel function of the 
third kind, first order, ^^lobai is the global driving force. 
The damping of the vortex motion rj is set to 1 (critical 
damping). Table summarizes typical values for the 
force and distance constants. 

To test the simulation for lattice distortion wc put a 
single defect within a toroidal area of 2 x 2 /im^ . The pin- 
ning force had a maximum value of -1 uN and a gaussian 
profile with a width of 25 nm. The driving force had a 
value of 10 pN leading to an average velocity of 9.8 pm/s. 
The vortex density was chosen to be ~ 242 /im~^ (near- 
est neighbor distance: ^ 69 nm ^ B — 0.53 T), similar 
to the vortex density within the data series considered 
in this paper. The simulations show similar velocity and 
distortion patterns as the measurements. Fig. 3] shows 
a map of the lattice distortions and the velocity around 
the defect. The distortion clearly reflects the attractive 
nature of the defect since all nearby vortices are deflected 
inward. Additional disturbances with decreasing magni- 
tude, showing roughly the lattice spacing and orienta- 
tion, can be found farther away. The velocity, however, 
mostly shows a repeating pattern of variation without a 
clear indication of the defect position. 

The different influence of a point defect on the dis- 
tortion and velocity is also found in distance and time 
dependent plots (Fig. [5] and IH]). The plots show traces of 
the average values as well as the envelope. The single de- 
fect influences the velocity almost uniformly throughout 



FIG. 6: Time dependence of the magnitude of the lattice dis- 
tortion and velocity. The plot shows the upper/lower bound 
as well as the average. The periodicity in the average velocity 
is largely due to the periodic boundary conditions. 



the simulated area leading to a modulation of the velocity 
in the time domain. In contrast, significant distortions 
arc only found within ~ 3 lattice Constance around the 
defect, while in the time domain the average distortion 
value is almost zero. Only the envelope indicates dis- 
tortions of up to 2 nm, but without showing a regular 
pattern. 

To find the decay rate of the distortion with distance, 
we fit an exponentially decaying periodic function: 

/(x)=A-e"^-sin(^7r-^-^^ (1) 

to the distortion in Fig. [S] using the least-x^-method. 
The parameters are given in the figure. The fit yields a 
period of 2 • w = 62 nm and a decay rate of Iq = 53 nm. 
The slight deviation of the periodicity is probably due to 
the attractive character of the defect. The decay rate is 
apparently not in an easy way related to the parameters 
of the simulation. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

Distortions of the vortex lattice were observed in all 
^ 50 time series we measured over the course of one 
year using 5 different NbSe2 samples at a temperature of 
4.2 K. The samples were about 5 mm in diameter and 
up to i < 0.5 mm in thickness. They were cleaved at 
a pressure p < 10~^ mbar before being introduced into 
the STM mounted inside the Helium cryostat. After the 
coarse approach the tip is roughly located in the center 
of the sample. Although we do not have optical access to 
the microscope, the distance of the tip from the center of 
the sample can be estimated. If we assume an expanding 
vortex lattice as the sole driving force, the distance to the 
expansion center r is given by r = ^^y^^ - a. a denotes 
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FIG. 7: Topography overlayed with the track patterns of the 
vortices passing through the area of three consecutive time 
series. The inset replicates the topography for clarity. The 
lines in the topography are single layer steps of NbSe2 (height: 
1.25 nm). 



the average lattice constant, v the average velocity. The 
value of r is typically below 1 mm. 

In this paper we focus on three consecutive time series 
taken within the same area of a sample (here: r = 0.15 
mm). The initial magnetic field of 0.5 T was stored in 
the superconducting magnet after introducing the sam- 
ple. The initial relaxation of the vortex lattice leads to an 
exponential decline of the average velocity with a decay 
constant of ~ 30 min. After a few hours a base velocity of 
in the order of pm/s remains due to a slow decay of of the 
magnetic field. The average vortex density shows a slight 
decline in magnetic field from 0.500 T to 0.497 T during 
the first time series with a finally value of 0.490 T after 17 
days. The average lattice constants and standard devia- 
tions are 69.2 ± 1.6 nm, 69.4 ± 1.5 nm and 69.8 ± 1.7 nm, 
respectively. In between each series, the tip was with- 
drawn from the sample surface in order to refill the liquid 
Helium dewar. Topographic images confirm that the lo- 
cation is identical to within one pixel {— 3.125 nm). Fig. 
[7] shows the topography as well as the combined vortex 
tracks. Several single steps are visible as well as a particle 
in the lower left corner and a depression running down 
the right side. Depressions like this one are regularly 
observed on layered materials. The tracks run through 
cluster points reproducing the vortex lattice. The points 
are connected by seemingly random tracks due to the 
complexity of the underlying defect pattern. However, 
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x-position [nm] 

FIG. 8: (Color online) Overview of the lattice distortion of 
three consecutive time series. The arrows represent the aver- 
age distortion within an 8 x 8 nm^ area. The three data series 
are plotted in different shades of gray. The circle marks an 
attractive site for the vortices. 



each time all vortices follow a particular track owing to 
the strong vortex-vortex interaction. The later is also re- 
sponsible for the high level of order in the vortex system. 
Although the vortex lattice should break up into domains 
we only observed six fold coordination. Calculating the 
distance dependent pair correlation functioriF^ shows a 
constant behavior with B{r) ~ 5.9 • 10~^. This is not 
entierly unexpected since we are only observing a small 
area of the vortex lattice of 1.6 x 0.4 ^m^ — including 
the lattice motion. 

It is still possible to pick out and study single defects 
by looking at the local distortion patterns. Fig. [8] shows 
an overview of the lattice distortion within the field of 
view. Each arrow represents the average distortion of 
an 8 X 8 nm^ area within the respective data series and 
points from the 'expected' to the 'corrected' position. 
The length of the arrows has been increased by a factor 
of 4 to improve visibility. Generally, the three data series 
show a good agreement, but there arc a few variations. 
To test the similarity we calculated the correlation of the 
local lattice distortion between data sets. The results 
are presented in table IIIII For the lattice distortion the 
values show a very high correlation (> 0.8) when com- 
paring series 1 and 2 and still a good agreement (> 0.5) 
for series 1 vs. 3 and 2 vs. 3. In comparison the veloc- 
ity shows little or no correlation. The differences might 
be explained by the slight increase in the lattice constant 
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0.53 
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0.39 


0.05 


0.35 


2 vs. 3 


1431 


0.61 


0.57 


0.58 


0.44 


0.18 


-0.01 



TABLE III: Correlation between data series for the 50 x 60 (~ 
8x8 nm^) raster. A'': Sample size, dx, dy: distortion in x and 
y, respectively, and {v±): velocity parallel (perpendicular) 
with respect to average travel direction. 
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FIG. 9: (Color online) Detailed view of the distortion pattern 
around subsurface defect. The arrows represent the average 
distortion within an 2 x 2 nm^ area. The three data series 
are plotted in different shades of gray. The arrows within the 
circle point to a common center. The length of the arrows 
has been magnified by a factor of 2. 



leading to a changing interaction with the defect pattern. 
It is noteworthy that the steps visible in the topographic 
images do not lead to any discernable effect on the vortex 
lattice. Previous studies using the Bitter technique^ii^ or 
scanning SQUID microscopj^ found clear effects due to 
steps. However, the studies where carried out at much 
lower magnetic fields (up to 3.6 mT) leading to a softer 
vortex lattice and for steps of at least 25 nm in height. 
In our case the step height represents less than 0.01 % 
of the sample thickness and is thus unlikely to have a 
measurable influence. 

A circle in Fig. [8] marks an area were the distortion 
is relatively high and is apparently pointing to a center 
position. Similar locations can be found, e.g. in the 
lower left corner. At the marked position, however, the 
topography shows no surface features such as steps or 
particles and was therefore chosen for closer examination. 
Fig. [5] shows the same location with a resolution of 2 x 
2 nm-. The average number of samples per arrow is 6.5. 
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FIG. 10: (Color online) Dependence of the magnitude of the 
velocity and the radial component of the lattice distortion on 
distance from presumed defect location. The plot shows the 
upper/lower bound as well as the average. 
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FIG. 11: Dependence of the magnitude of the velocity and 
lattice distortion on time. The plot shows the upper/lower 
bound as well as the average. 

The length of the arrows has been magnified by a factor 
of 2. The arrows within the circle point to a common 
center where we suspect a subsurface defect. We also 
note, that the defect hes within a cluster point and has a 
focusing effect as the vortices pass by. The motion occurs 
essentially from the lower right to the upper left corner. 
Still, there are variations in the vortex tracks, even in 
close vicinity of the defect. Therefore, while this pinning 
center appears to have an infiuence on the vortex motion, 
other pinning sites have to be involved to form the track 



and velocity patterns observed. 

Finally, Fig. [TU] and [TT] show the dependence of the 
lattice distortion and velocity on the distance from the 
center found in Fig. IH] and on time, respectively. The 
curves show a similar behavior to that found in the sim- 
ulations (cf. Fig. [5] and [5]) . The distance dependence 
of the distortion shows a clear first maximum and oscil- 
lations, while the velocity is featureless. The roles are 
reversed when analyzing the time dependence. The ve- 
locity shows several peaks, stemming from the collective 
interaction of the vortex lattice with the defects, whereas 
the magnitude of the average distortion is essentially flat. 
Fitting equation [1] to the distance dependence of the lat- 
tice distortion yields a period of 2 • m; = 78.8 nm and 
a decay rate of to = 291 nm. Given the large width of 
the envelope after the first maximum and the poor cor- 
relation coefficient of = 0.46, one cannot give much 
weight to these numbers. 



V. CONCLUSION 

We discovered local lattice distortions in slow moving 
vortex lattices in NbSe2. The distortions were found in 
2D simulations as well as in measurements. It is a lo- 
cal effect generated by the interaction of vortices with 
individual pinning sites. Thus, observing the distortion 
patterns allows to identify the location. This is the first 
step for a detailed study of the vortex-defect interaction. 
In contrast, variations in the velocity are a collective phe- 
nomenon generated through the collective interaction of 
the vortices with a larger group of defects outside of our 
observation area, facilitated by the strong vortex-vortex 
interaction. This becomes evident by examining the dis- 
tance and time dependence of the distortion and veloc- 
ity. While the distribution of defects governs the vortex 
tracks and the patterns in the velocity, the elasticity of 
the vortex lattice still permits individual vortices to be 
influenced by individual defects. 
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